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Abstract 

The motion of a conducting electron in a quantum dot with one or several 
dislocations in the underlying crystal lattice is considered in the continuum 
picture, where dislocations are represented by torsion of space. The possi- 
ble effects of torsion are investigated on the levels of classical motion, non- 
relativistic quantum motion, and spin-torsion coupling terms derivable in the 
non-relativistic limit of generalizations of the Dirac equation in a space with 
torsion. Finally, phenomenological spin-torsion couplings analogous to Pauli 
terms are considered in the non-relativistic equations. Different prescriptions 
of classical and non-relativistic quantum motion in a space with torsion are 
shown to give effects that should in principle be observable. Semi-classical 
arguments are presented to show that torsion is not relevant for the classical 
motion of the center of a wave packet. The correct semi-classical limit can 
instead be described as classical trajectories in a Hamiltonian given by the 
band energy. In the special case of a spherically symmetric band this mo- 
tion reduces to straight lines, independently of local crystal orientations. By 
dimensional analysis the coupling constants of the possible spin-torsion in- 
teractions are postulated to be proportional to a combination of the effective 
mass of the electron, m e g, the lattice constant, a, and fi. The level-splitting is 
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then very small with transition frequencies on the order of 1 kHz or smaller. 
PACS numbers: 61.72Bb, 73.20.Dx, 85.30.Vw 

I. INTRODUCTION 

Quantum dots are small essentially two-dimensional conducting domains, connected to 
external leads by tunnel barriers such that the number of conducting electrons is quantised. 
The elastic mean free path and the dephasing length may, under appropriate circumstances, 
be much larger than lateral size of the dots. The motion of an electron inside may therefore 
to a first approximation be considered as ballistic. These structures may hence in many 
respects be regarded as physical realizations of the quantum mechanics text-book examples 
of motion in two-dimensional potential wells. For recent reviews and papers, see 

The purpose of this Letter is to explore a suggestion made in a book and in a series of 
papers by Kleinert and co-workers P~PTJ|, which I will state compactly as follows: i) the 
motion of a quantum particle in a space with torsion has several distinct features compared 
to free motion in flat space; ii) a crystal with defects may in the continuum picture be 
described as a space with torsion; iii) the motion of an electron around a defect may thus 
in some circumstances be modelled by the motion of a quantum particle in a space with 
torsion. 

A continuum description of a crystal is valid at distances much larger than the lattice 
spacing. The low- lying electronic states in a quantum dot extend across the dot, with 
typical spatial scale L, while a localized defect extends over a few lattice spacings a. The 
scale separation L/ais thus in the range of 100 — 10000, and it is conceivable that an effective 
continuum description might be appropriate. 

The outline of this paper is as follows. In section [Tl] I review the literature on motion in 
spaces with torsion. The Dirac equation coupled to torsion allows for one vector and one 
axial vector coupling, both of which give rise to Zeeman-like couplings of spin to an axial 
vector in the non-relativistic limit. We also consider more phenomenologically motivated 
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spin-torsion couplings which yield similar terms. In section |T| I review elementary aspects of 
the differential geometry of physics of defects and in section [TV] I discuss on the semi-classical 
level the effects of crystal torsion on elecron motion. In section [V] I try to evaluate possible 
torsion-induced effects in quantum dots. In section [VT] I summarize the results. They are 
are negative as concerns the proposal by Kleinert that autoparallels should play a privileged 
role. On the semi-classical level it is possible to observe motion along geodesies, but such 
motion is insensitive to torsion. Motion along autoparallels and a related non-relativistic 
quantum mechanics could also probably be disproved by experiments. The spin-torsion 
coupling terms are not theoretically unsound, but leads to very weak effects that they would 
be difficult to observe. 



II. CLASSICAL AND QUANTUM MOTION IN A SPACE WITH TORSION 

In this section I use the Einstein conventions of summing over repeated indices and 
raising and lowering indices by the metric tensor and its inverse g % K For conventions 



pertaining to the affine connection, torsion and contorsion I follow Hehl et al [12|, see also 



Schrodinger |13| and Schouten [|14 



A manifold is said to carry metric and affine structure if one can compute the length of 
vectors in the tangent space by 

\A\ 2 = g^A* (1) 

and parallel transport of a vector along an infinitessimal distance dx by 

A i _> A* - F kj i A j dx k (2) 

The metric and the affine connection T are connected by the compatibility condition, that 
the scalar product of two arbitrary vectors is invariant under parallel transport: 

9ij,k — Qij^ki 1 — 9u^kj 1 = (3) 
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which may also regarded as the statement that g is covariantly constant. The first term in 
(0) stands for the partial derivative dgij/dx k . 
The standard form of the solution of ([$]) is 

da^kj 1 = 9u{^j } ~ Kkji ^kj^ = ^ mj ' k ^ mk ' j ~ 9kj,m) (4) 

where we recognize the Cristoffel symbol symmetric under the interchange of j and k. 
The second term in (H), K^, is known as contorsion and must be antisymmetric under the 
interchange of the last two indices of j and i. 

Let us consider a small area spanned by two vectors dxm and dxw. If we transport dx(o) 
along c?X(i) it will change into (in component form) cfo| 2 ) — Vj t j l dx 3 , 2 -,dx k x y If we move first 
along dxm and then along dxm we will hence end up with a displacement (in component 
form) 

dP^2 ~ dx ~t~ dx — ^kj ^dx^2\dx k ^ (5) 

If we were on the other hand to make the displacements in the opposite order we would end 
up at 

dl^2i ~ dx '(2) ~t~ dx^^ — r f^j dx^^dx^ (6) 

The difference between these displacements is a vector, which measures by how much the 
circuit fails to close if direction vectors are parallel transported around the perimeter of the 
area spanned by dxm and dxr 2 ) 

dB = dP^2 d^2\ ~~ ^- kj (dx^^dx dx^^dx^^j (^0 

The parenthesis in ([j]) is the area element dA\2^ k spanned by dx^i) and dx^)- The part 
of r antisymmetric in interchange of j and k is called torsion, and has hence the following 
relation to dB (in component form): 

dB 1 = —S kj l dA 12 kj S kj 1 = T[kj] 1 (8) 

Torsion is therefore a third-order tensor, and it may be connected to contorsion by 
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K-ijk Sijk -\- Sjki Sfcij (9) 

from which it is seen that the required antisymmetricity of K (in the last two indices) follows 
from the antisymmetricity of S (in the first two indices). As a final consistency check one 
may antisymmetrize V using the decomposition of (^): 

i 1 / 

T [kj] 1 = { ^.i } - K [kj] 1 = - - (K kj 1 - K jk 

= 2 ( Sk i 1 ~ Sjk l ) = S k j 1 (10) 

The symmetric part of T k j 1 under interchange of k and j thus contains both the Cristof- 
fel symbol and a symmetrized combination of the antisymmetric part, which explains the 
distinction between torsion and contorsion. 
The geodesies on a manifold are given by 

d oc r % dx^ doc / i \ / \ 

H^ + { ]k } 77!7 = (geod<islcs) 

The intuitive concept of free motion seems however to imply that velocity vectors change 
according to the law of parallel transport, and such curves are called autoparallels: 

d*^ cc^ dx^ dx^ 

—— + Y kj l — — = (autoparallels) (12) 

CLT 1 dr GST 

The interesting observation stressed by Kleinert is then that in a space with torsion geodesies 
and autoparallels do not coincide. It is not obvious which of the two is the most natural 
extrapolation from ordinary classical mechanics. For different points of view, see, on the side 



of autoparallels, Kleinert H-p], and on the the side of geodesies Hehl et al [12], Audretsch 



and the recent papers of Lammerzahl [24] and Barros e Sa ||15|| . 



In any case, from the viewpoint of Hamiltonian dynamics autoparallels have peculiar 
features. One of the simplest examples is motion on a two-dimensional surface with constant 
diagonal metric and constant torsion. The torsion tensor has then only two independent 
non-zero components, namely S'121 and 6*122, which together specify a direction in the plane, 
b = (5*121, 5*122) • It is straightforward to integrate equations (|T2"D: they here describe motions 



at constant speed, but where the direction of velocity changes so as to be orthogonal to b. 
In other words, kinetic energy in the ordinary sense is conserved, but momentum and phase 
space volume is not. In these respects motion along autoparallels in a space with torsion 
is similar to a mechanical system with non-holonomic constraints |l6|Jl7[| (for analytic and 
numerical investigations of a clarifying concrete example, see |l8|-p0[|), an analogy also used 
extensively by Kleinert. 

The difference between autoparallels and geodesies carry over to the quantum mechanics 
of a non-relativistic scalar particle. The gradient of a scalar is a vector, and the covariant 
derivative of this vector involves the connection, so that the contraction of two covariant 
derivatives acting on a scalar is 

g ij DiD^ = gVDidrf = (g l] d l d J - g^Ty l d t ) $ (13) 

The usual generalization of the kinetic term in the Schrodinger equation involves on the 
other hand the Laplace-Beltrami operator 

^= (d^^gdj) i, = - 9 l3 {\ 2 }dij i> (14) 

which is self-adjoint with respect to the standard scalar product < r}\ij) >= J ^fgrf^). 
The difference between (|l3~) and (14) is a gradient 



K^d^ = -2S l d^ S^Sn* (15) 

which, as remarked by Kleinert, is not in general self-adjoint, at least not with respect to 
the same scalar product. Quantum dynamics governed by (0) is therefore not necessarily 
unitary. This somewhat disturbing property can be compared with the fact that phase space 
volume in the usual sense is not conserved by classical motion along autoparallels. 

The problem of geodesies vs. autoparallels can finally also be considered, albeit a little 
indirectly, for relativistic spin-- particles. To define spinors in General Relativity one needs 
a system of local intertial frames each defined in the neighbourhood of some space-time 
point X. The transformation matrix between £j£ and a chosen, in general non-inertial local 
coordinate system at X is called a vierbien 



The first index of the vierbein stands for a direction in the local inertial frame, and is hence 
lowered and raised with the Minkowski metric rj. The second index stands on the other hand 
for a direction in the tangent space at X, in a basis determined by the coordinate system 
x^, and is therefore raised and lowered with the metric tensor g. 

Parallel transport and covariant differentiation of a spinor in direction a of the local 
inertial frame are determined by the vierbiens as 

= + W r^^v^W ( 17 ) 

where cr^ 7 are the basis elements of infinitessimal Lorentz transformations in the spinor 
representation, and the semicolon in the last factor stands for covariant differentiation of 
the vierbien according to 

V 7 v., li = V rV)lt -{* v }V yK (18) 
The Dirac equation in a space with curvature is then 

ih^ a V a ip + mclip = (19) 
For a classical discussion of these matters, see E1J. For a thoroughly modern discussion 



from the mathematical point of view, see Wl 



When we allow also torsion of space we enter into less well-chartered territory. One 
possible way to proceed is to take as building blocks the vierbeins directly, allowing for V a ^ v 
not necessarily being equal to V a „ )|t . The minimal formal change in the Dirac equation is 
then to perform the covariant differentiation in (|18D using the full affine connection instead 
of the Cristoffel symbol. However, we have now also torsion as an independent tensor. By 
contractions we can form a vector S = S l m m and an axial vector S = e no S mno . The 
simplest relativistically covariant additional terms in (|T9|) are hence 



(iC%S al a + DhS al5l a ) $ (20) 
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where C and D are dimension- less numerical constants. In fact, doing the covariant differ- 



entiation in (Tl8| ) with the affine connection will not bring in any more terms except those 
in (^). The Dirac equation coupled to torsion by the simplest vector and axial vector 
couplings and covariant differentiation is thus 

ifi^VJ) + ((C - l)ihS a ~? a + (D+ h)hS al5l A ip + mclt/j = (21) 

where the covariant derivative T> a is defined by (|17|) and CPL8|) . If we also put the constraint 
that (|2~TD should be derivable by variation from a real action, we are led to the specific choice 
C — 1, D — 0, which is often in the literature referred to as the Dirac equation minimally 
coupled to torsion [ ^ , |T2|f2~3| , f24}1 . 



I will in the following neglect time components of the space-time torsion vector Si. The 
only non- vanishing component of the axial vector Si is then So, which, if we look at it in 



3D, transforms as a pseudo-scalar. The axial vector coupling in fl2~T|) has then the following 
non-relativistic limit [[24[] 

i(D + h-S (a-d)i; (22) 

where a is the vector of Pauli matrices and the derivative acts on the wave function to the 
right. 



The vector coupling in (£L|) is on the other hand completely analogous to the coupling 
to an external electromagnetic field. It will thus, in the non-relativistic limit, give rise to an 
effect like the coupling of spin and magnetic field 

^-(Vx5)-# (23) 
2 m 

where S stands for the vector of spatial components of the four-dimensional torsion vector. 

In summing up this section we see that the equations for autoparallels flT2| ) and the 
putatitive generalization of the Schrodinger operator as the contraction of two covariant 
derivatives formed with the affine connection ( pT3|) are the odd ones out. Neither arises as 
the limit of the Dirac equation minimally coupled to torsion, which only gives a spin-torsion 
coupling reminiscent of a Pauli term, the standard Laplace-Beltrami operator and geodesies. 
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III. TORSION: CONTINUUM REPRESENTATION OF DEFECTS 



The description of defects in crystal lattices in the continuum picture with tools from 
differential geometry has a distinguished history, described in the 1980 Les Houches lectures 



of Kroner [26], which also give a good introduction. 

The basic idea is quite simple. A crystal is supposed to carry only isolated defects, such 
that one has around most points a locally perfect lattice. The lattice locally around a point, 
if continued without defects or deformations indefinitely, may be considered as the tangent 
space of the crystal at this point. The local lattice directions provide a basis for this space. 
A vector in the tangent space can then be identified with moving a certain number of lattice 
units in each direction. Parallel transport of a vector from lattice point P to lattice point 
Q means that we identify the vector of, say, rii steps in crystal direction i at P with the 
vector of an equal number of steps in the same lattice direction at Q. Since the local lattice 
directions may change from point to point in the crystal, a vector can change under parallel 
transport if measured in an external frame of reference. 

We can then consider the following process. Take two lattice vectors n\ and 112 in 
directions %\ and %i- Transport the pair first along m, then along 712, —n\ and — ri2. If the 
crystal is perfect the circuit closes and we are back at the point where we started. 

If, however, the circuit has circled a dislocation line, we are not back at the point of 
departure. Dislocations lines in three-dimensional crystals are characterized by a vector t 
tangential to the line, and a vector b describing the mismatch if we circle the dislocation 
line in the positive sense determined by t. 

That is, if we introduce a vector dAj normal to and equal in length to the area spanned 
by rii and 712, then the mismatch when going around the circuit is linearly related to dA by 

db l = a ji dAj (24) 

— * 

where a is the density of dislocation lines with t in the j direction and b in the i direction 
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The procedure described here is of course identical to that used to define the torsion 
tensor in section |IT] so one may introduce a crystal torsion field as 

Skj 1 = ~ -tkjmOt" 11 (25) 

where e is the totally antisymmetric Levi-Civita tensor and summation of repeated indices 
is understood. 

In three-dimensional simple crystals with one atom per Bravais cell there are two quali- 
tative types of dislocations: screw dislocations where b is parallel to t; and edge dislocations 



where b is normal to t |2lf-|28|j. In the first case the dislocation tensor a has only diagonal 
elements, while in the second case it has only off-diagonal elements. 

The totally antisymmetric trace of the crystal torsion field vanishes for edge dislocations 
since 



e k3t S k3l = -a\ (26) 

On the other hand, for screw dislocations the contracted torsion vector vanishes since 

Si = Sii 1 = —-eu m a ml (27) 

Si can hence be looked upon as a vector dual to the antisymmetric part of the dislocation 
density tensor. For edge dislocations it is normal to both b and t, and equal to half of b in 
length. 



IV. GEODESICS VS. AUTOPARALLELS IN CRYSTAL SPACE 

A crystal space with dislocations but without interstitial defects carries torsion but not 
curvature |26| . Geodesies in this space means simply straight lines in the frame of reference 
of an external observer. Autoparallels means on the other hand motion which is always in 
the same direction with respect to the local crystal directions. 

The following discussion will be on the level of a semiclassical model. The basic idea is 
to take a wave packet, a superposition of Bloch states, and assume that the crystal changes 
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on length scales much larger than the spread of the wave packet. To construct the Bloch 
states I can therefore take the crystal lattice close to a given point where the wave packet 
is centered and extend it indefinitely in all direction without deformations or defects. This 
corresponds to the tangent lattice in the sense of section [TTT| . If I go far enough from the 
point the real lattice and the tangent lattice will differ, but at those distance the amplitude 
of the wave packet is assumed vanishingly small, so this difference will be ignored. 

More quantitatively, the wave packet is built from Bloch states with wave vectors in a 
domain of size Ak around k. We assume that Ak is small compared to the dimensions of 
the Brillouin zone, and that 1/Ak is small compared to the scale L on which the lattice 
changes orientation appreciably. In real space the wave packet is then located in a domain 
of size 1/Ak around a center, which is denoted r. 

In the absence of external electric and magnetic fields the equations of motion for the 



center of the wave packet in an undistorted crystal are [27 



r = v„(k) = ^^ hk = (28) 

where hh is the crystal momentum and S n (k) is the energy of the state with wave vector k 
in the n'th band. Equation Q2lf ) describes motion in a straight line. This picture is valid as 
long as dispersion effects are not important, 

* « *dis P - sjsJ {Ak)2 ( 29 ) 

To see the effects of the changing crystal before dispersion sets in too strongly we must 
demand that L is much less than the distance the wave packet traverses during time tdisp> 
that is 

Y d£„ (k) ^ 

— « L « L disp L disp = a 2 ff (k) ( 30 ) 

ak 2 v ; 

As we will see later, for realistic values of Ak the bounds in ([3(]) are a little tight, but for 
the present discussion it is sufficient that there are in principle some scales of time where the 
wave packet is still localized on a length scale much less than L, but has moved a distance 
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much larger than L. In this intermediate regime one can thus pose the problem if the wave 
packet would follow geodesies or autoparallels or some other curves in crystal space. 

Let the motion of the center of the wave packet be parametrized by (x*(t),k(t)), both 
from now on given in a coordinate system fixed in space. The changing orientation of the 
crystal can then be described by postulating that the band energy £ n (k, x) is a slowly varying 
function of x. As the wave packet moves in the crystal from position x*(0) at to position 
x*(t) at t it acquires a phase of e % ' ( »' ( " where 

S[x*(t),k(t)} = [\hk(t)-x*(t)-£ n (k(t),x*(t))} (31) 
J o 

The actual path is determined by the condition that S should be stationary under variations. 
It is obvious that if we introduce the momentum p = Tik and the Hamiltonian function 
H(p,x) = £ n (jr,x) the center will follow the classical trajectories of H. 

Let us assume for simplicity that the band energy has the structure £ n (k) = \% 2 m~^k l ki , 
where m^ is the effective mass tensor, which has principal axes , . . . , and effective 
masses in those directions mi, ... , m^. The changing orientation of the crystal is effected by 
letting the principal axes depend on x. The classical Hamiltonian is then 

E 2^r(P • ^(*)) 2 = E ^(*)PV 9ik(x) = E ^-nf{x)nf{x) (32) 

i 1 jk i 1 

The wave packet will hence follow the geodesies with respect to the metric induced by 
the crystal orientations and the band structure. 

In the special case of a spherically symmetric band structure the effective mass tensor is 
proportional to the identity. The Hamiltonian is then ^m^p 2 , and the dependence of the 
crystal orientations drop out. In the last idealised case we would thus predict motion along 
straight lines which are geodesies and not autoparallels in crystal space. 

To end the discussion we must also take into account that lattice directions in real crystals 
are changed by the presence of dislocations. These are local scatterers, and we have to check 
if the coherence of the wave packet can be maintained over such distances that the lattice 
orientations change appreciably. Suppose that the strength of the perturbation is U and 
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that it has support over a typical distance of the lattice spacing a. The wave packet passes 
over the defect during a time T scatter of about — Je n (k) > while the natural time scale of the 
action of the scattering potential on the wave packet is h/U. If we assume scattering to any 
state on the energy shell with equal probability, we have 



^scatter on shell ~ [dAk] (ak) (^^J ~^~^J ( 33 ) 

Equation ( [33]) is an overestimate, but probably not very much so. A derivation of (^3|) using 
standard first-order perturbation theory is given below in appendix 

Consider now an array of dislocation lines with surface density n and a wave packet 
moving in a plane perpendicular to the lines. As it traverses a distance L it will in general 
encounter L-^n defects, and the total probability of being scattered by any of them is 

Scatter along L ~ (l^j-n) (aAk) d W (jf) (^) 

We want this probability to be much less than one. A surface dislocation density of n leads 
to a crystal torsion field of strength na, since each dislocation contributes a Burgers' vector 
of length a. The crystal directions can change appreciably over a length L if the line integral 
of torsion is of order one, that is naL ~ 1. The two estimates are compatible if 

(^i"^)'^) (f) <<x (35) 

We are interested in the case d = 2. Assuming for simplicity £ n (k) ~ |m off k 2 we have 

(ok)3 (m) 2<<1 (d=2) (36) 

This bound can evidently only be satisfied, and then only for sufficiently large crystal mo- 
menta k, if the effective scattering potential U is much less than the the largest kinetic 
energies in the band. Condition fl3ll| ) is a serious point. If U is not sufficiently small then 
scattering is the leading effect and considerations pertaining to fine points of the continuum 
picture, which we pursue here, are simply irrelevant to electron motion in quantum dots. 
We remark in contrast that in the 3D situtation (BBl) becomes 
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u 



) 



2 



(aAk)(ak) 4 



S n (k) 



« 1 



(d = 3) 



(37) 



which can always be satisfied if Ak is small enough. 



V. POSSIBLE TORSION EFFECTS IN QUANTUM DOTS 



I will now consider the following possible types of effects: i) classical motion; ii) 
non-relativistic quantum motion; iii) spin-torsion coupling terms derivable from the non- 
relativistic limit of the Dirac equation minimally coupled to torsion and iv) phenomenolog- 
ical spin-torsion couplings appropriate for the situation of quantum dots. 

When classical or semi-classical descriptions are valid the electronic motion in the dot 
is effectively two-dimensional. Let the two directions of the dot be 1 and 2 and the third 
vertical direction 3. For points i) and ii) we only need to consider the possible effects of 
elements Swi and S122 of the crystal torsion matrix, i.e. only edge dislocations, while for 
iii) and iv) we also need to consider screw dislocations. 

It is convenient to express torsion in lattice unit a. Kroner discusses a crystal with one 
edge dislocation per every ten-by-ten atoms, all Burgers vectors of all dislocations oriented 
the same way. In this situation crystal torsion would be 0.01/a, and one would have to move 
a distance 100a to see to change the lattice directions change appreciably. It is not clear if 
this example is realistic, the surface defect density being about 10 17 m -2 . As the standard 
set-up I will thus consider one a dislocation per every hundred-by-hundred atoms. Torsion 
would then be 10 _4 /a, and one would have to move a distance of 10 4 a to change the crystal 
orientations. Crystals with lower dislocation densities could surely be fabricated, but then 
we would need macroscopic quantum dots to see the possible effects of crystal torsion. 

With quantum dots of micron size the bounds on dispersion from section [IV] can be 
satisfied. We would see the lattice orientations change when we move from one side of the 
dot to the other, i.e. L ~ 10 4 a. If we choose the spread of the wave packet Ak to be about 
10~ 3 /a. then 1/Ak << L « Ldisp with an order of magnitude at the lower inequality 
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and two orders of magnitude at the upper inequality. The semi-classical argument predicts 
that wave packets follow classical trajectories in a Hamiltonian given by the band energy as 
functions of crystal momentum. This structure stems from solving for the Bloch states in the 
crystal potential, an information which is not contained in the continuum description of the 
crystal as a metric space with curvature and torsion. In the general case neither geodesies nor 
autoparallels in crystal space therefore have any particular relevance to the problem. In the 
special case of a spherically symmetric band structure we do however recover motion along 
straight lines, i.e. geodesies in crystal space. Motion along autoparallels would on the other 
hand be on curved paths as in the example discussed in section ||. If the approximation of 
a spherically symmetric band structure is a sufficiently good one, and if the local scattering 
potentials are weak enough that the wave packet does not loose coherence, then the difference 
between geodesies and autoparallels could in principle be observable. 

The scale separation is even more favorable for a description in terms of a non-relativistic 
quantum particle. Consider the ground state or a low-lying excited state in the example 
discussed before. If the prescription flIB|) were correct, they would be solutions to the 



eigenvalue equation 

,2 *2 



h 1 



- V 2 +yu(x- xA + S l d t U = Ei/j 



(38) 



2m eff 

where U is the scattering potential from dislocations at points Xi and S l is the crystal torsion 
vector, which in the example under discussion would be about 10 _4 /a. The derivative d\ 
acting on a low-lying state is also about 10 _4 /a if the lateral size of the dot is in the micron 
range. Hence the perturbation S di is of order 10~ 8 /a 2 , which is comparable to the energy 
gap between the ground state and the low- lying states of (p8|). As discussed in above in 
section || the perturbation breakes hermiticity, and hence leads, if large enough, to complex 
eigenvalues, i.e. states exponentially decaying or growing in time. On a conceptual level the 
problem is primarly the exponentially growing modes. One might imagine that the defects 
could perhaps excite certain electronic states in the dot, and that for some transient time 
the growth rates of those states could be decribed by the eigenvalues of fl38"l) . This picture 
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does not look physically very plausible, and could probably be ruled out by experiments, at 
least on a qualitative level. 

Let us now turn to possible spin-torsion couplings. The most straight-forward, but also 
certainly the smallest, are the non-relativistic limits of the vector and axial vector couplings 
in the Dirac equation, ( P3| ) and (^). In the first case we have the curl of the torsion vector. 
Let us assume that the torsion varies in a controlled manner in the plane of the dot. We 
only have to consider edge dislocations, since the vector S vanishes for screw dislocation. 
When the density of edge dislocations varies normally to the vector S we would have an 
effective coupling analogous to the Zeeman-effect, with an effective magnetic field pointing 
in the vertical direction (out of the plane of the dot). In the axial vector case we only have 
to consider screw dislocations, S being the sum of densities of screw dislocations in all 
directions. 

We can also discuss more phenomenological spin-torsion coupling terms. In the physical 
situation of a quantum dot the vertical direction is determined by the gradient of dopant 
concentration. We can therefore form a combination of torsion and the vertical direction 
vector which transforms as an axial vector, namely S l = ej ^ k n\Sjk 1 ■ This is just the vector 
b which results from circling a dislocation in the 12-plane, i.e. S could have elements 
both in and out of the plane. We then consider a term like kS ■ a. Torsion has physical 
dimension of inverse length, hence k should have dimension mass • (length) 3 (time) -2 . The 
only combination of Planck's constant H, the lattice constant a and the effective electron 
mass m e g, which has this dimension is . This coupling can be compared dimensionally 
with ( p3[) and (p2|), which have the bare electron mass instead of m e s and a space derivative 
instead of 1/a. 

The energy splitting of two states with spin respectively parallel and antiparallel to S 
will then be will then be k times the average value of S in lattice units. With the standard 
case of one dislocation per hundred-per-hundred atoms this gives 

Aversion « ^ ^ 
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The frequency of the transition is hence 

torsion ~ 10~ 4 ^-^ ~ 100Hz (40) 

m efi a 2 

where I have for simplicity estimated m e s with the electron rest mass and a to be 1 A. The 
two terms derivable from the Dirac equation discussed above are even smaller, approximately 
by a factor a/L. 

VI. CONCLUSIONS 

I have in this paper taken up a suggestion of Kleinert that the motion of a classical of 
a quantum particle in a space with torsion could be relevant to describe some properties of 
crystal with defects. I have focused on the example of electron motion in quantum dots, 
since these systems most nearly realise the textbook example of electron motion in potential 
wells. 

In a space with torsion there are two different priviliged classical motions, geodesies and 
autoparallels. A semi-classical analysis carried out in section [TV] indicates that at some 
scales of space and time geodesic motion could be observed, if the band energy structure 
is spherically symmetric, and the local scattering potential from dislocations sufficiently 
weak. The last condition is very stringent in 2D geometries such as quantum dots, and local 
scattering is in fact likely to be the leading effect. A continuum description of the motion 
of an electron in a crystal with dislocations is therefore probably not very useful in 2D. 
Assuming nevertheless, for the sake of the argument, that the local scattering potential is 
weak, the conditions for the distinction between geodesies and autoparallels could probably 
be experimentally realised. These issues are discussed in section [V|. If the band energy 
structure is not spherically symmetric the problem is more complicated and neither geodesies 
nor autoparallels would follow from the semi-classical argument. 

The distinction between geodesies and autoparallels carry over to the quantum mechan- 
ics of a non-relativistic particle, i.e. to the proper generalization of the Schrodinger equation 
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to a space with torsion. The analogy with geodesies would be the Laplace-Beltrami opera- 
tor, while the analogy with autoparallels would be another differential operator introduced 
by Kleinert. The difference between the two operators can be relatively large, and hence 
observable. Since it breakes hermiticity. it leads to modes exponentially growing or decaying 
in time, an effect which should be falsifiable by experiments. 

The spin-torsion coupling terms, both the ones derivable from the Dirac equation and 
more phenomenological ones, are all very small and give transitions with frequencies less 
than 1 kHz, or even less than that. 

Hence we conclude that whatever the motion of a classical or quantum particle in a 
space with torsion should be, the issue has little bearing on electron motion in quantum 
dots. The semi-classical motion of wave packets can, in favorable cases, be arranged to 
be along geodesies, i.e. along straight lines, but autoparallels are never relevant. Bound 
electronic states in a quantum dot may to some accuracy be described as the eigenstates of 
the Laplacean with Dirichlet boundary conditions, but the alternative introduced by Kleinert 
has eigenstates that are exponentially growing or decaying in time, and this contradicts 
the very notion of stationary bound states in the dot. The spin-coupling terms are not 
in contradiction with theoretical arguments, nor with concievable experiments. But their 
effects are so weak it would be very difficult to observe them. 
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APPENDIX A: COMPUTATION OF THE WAVE PACKET-DEFECT 

CROSS-SECTION 

The wave packet interacts with the defect during a time of about 
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It is therefore convenient to regard the defect as a time- dependent perturbation, switched 
on before — T scat ter and switched off after T scat ter- For definitiveness take V(x,t) = 
V(x) exp(— 2^(£/T scatter ) 2 ) with a some number which will eventually be taken to infinity. 

Let us take a lattice of large finite size A with the isolated defect centered at the origin. 
The values of the wave vector k then come in integer multiples of 2tt/A. The time-dependent 
Bloch states are denoted ty^\(x), and are normalized by J Ad |^ / i°i(a ; )| 2 = 1- The wave packet 
formed by Bloch waves from the n'th band is given by J2k ^n\ a k°\ where the amplitudes 
are normalized by J2k \ a k°^\ 2 = 1- 

When the wave packet has passed over the defect, i.e., after the perturbation has been 
switched off, we have to first order in the perturbation the out state 

¥(out) = £ *2 (4 0) + £ U kl af^ U kl = -\ Vu{ty^ (A2) 

where Vki(t) is the time-dependent matrix element between the time-independent Bloch 
states k and I. The time dependence of the Bloch states have been written out explicitly 
in the matrix element in ( |A2|) . The difference between ^(out) and the unperturbed wave 
packet at the same time is 5^ (out) = J2ki ^i°l^fcz a z°' ) - It is orthogonal to the unperturbed 
state. The absolute square of S^f hence represents the probability to be in an out state 
orthogonal to the unperturbed wave packet. 

I disregard scattering out of the band. Then, for pairs of states such that \E k — E{\ » 
Tp — the transition probability amplitude Uki will be very small. Inserting the definition 

-'scatter 

of T sca tter this means that only wave vectors in a thin shell of thickness about Ak around 
the surface given by = E^ are scattered. The matrix element is 

V kl (t) =< k\V\l > exp(-^(t/T scatt er) 2 ) (A3) 

where < k\V\l >= J (^^ k )* (x)V (x)^](x) dx. If both wave vector are within Ak of k then, 
because 1/Ak is assumed much larger than a, 
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< k\V\l >= (~\* U \k-l\~Ak (A4) 

where U is some measure of the strength of the potential. In general, if \l > is in the wave 
packet we can always for the same reason substitute < k\V\l > with < A;|V|k >. We have 
now 

% ( 1 / „ ,E h -E^^ 



U kl = -- < k\V(x)\l > exp ^T scattcr (^ 7 -^)J J v /2 7 ra 2 T s 2 catter (A5) 
The total probability of scattering is 

i^i 2 = E^M4 0) r«S 0) (A6) 

Let the amplitudes be given as = fo(k), where the smooth function /o(&) is 

normalized by J Ad |/o(&)| 2 = 1- For defmitiveness take fo(k) to be a Gaussian centered at k 
with width Ak: 

1 ( k - k\ 2 \ / / \ t ^2\~ d / 4 



Uk) I ^[-^-j J (nW) (A7) 

We now wish to compute J2i U n iaf^ where n is anywhere on the shell of scattered wave 
vectors. We have to separate components of wave vector / parallel and orthogonal to 9£ q^ ■ 
By expanding 

E n -E l = (E n - E k ) - ^^-(h ~ k,|) (A8) 



and using the definition of T scat ter from ( |A1| ) we have 



E = ~< n\V(x)\k > p^T^ ter (— j (^(A^f (27r(Afc) 2 



A\ d/2 / . A ,,o\-<i/4/ , AIN o N \<*/2 



1 / 1 CT 2 {T S c&tter(E n — Ek) X 



1 1 scatter v-^n -"K/ i i /a 
eX P -0T2-TT i ( A9 ) 



V tr 2 + 1 ^ ^ 2 a 2 + 1 V ft 
The contribution to ( |A6|) from a given wave vector n can therefore be written 

< nlVlk > |T sra f tnr \ 2 a 2 



E^^H 'V ) ^ rT (AA;) d (2 7 r)^ +1 2^ 



ft / (T 2 + l 
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When we sum over n we can take out a value n such that E n = E^, and then integrate 



parallel to That will give a term J jy^|w>j • Incorporating the partial derivative 

V scatter Sn > 

in a delta function we can hence write 



= /_^(AA;) d + 1 (27r) (i / 2+2 2 ci / 2 - 1 (Te^l) \?E 
V (J + 1 V ft / 9k 



2^ 



f8(E n -E k )\ < n\V\k> \ 2 dn (All) 



Here we can take the limit of a going to infinity and we have thus removed the spurious 
time-dependence of the perturbation. 

By order-of- magnitude estimate we take | < n|V|k > | = U (j^) for all states n on the 
energy shell given by E n = E^. We can then take the matrix element outside the integral 
/ 5(E n — Ek) dn, which is the density of states at energy E^. If we further assume that 
E n = ^m c gn 2 , then the density of states is equal to ^fjr-, where Qd is the area of the unit 
sphere in d dimensions. We can finally take out from ( [All] ) the combination Ak\^- |^A- 
which can be rewritten more simply as and this gives 

If we assume instead that matrix element is only non-zero between states in the wave 
packet, then we can estimate the integral in (|A11|) by ji n _ k i<Afc &{E n — E k ) dn which is about 
(Ak)^ 1 The probability of being scattered from one wave packet to another is hence 



"Wain r u-il.liiii A A ~ (//A//)"" ( "" : "' r ) (A 12) 



which can also be derived by straight-forward dimensional analysis, since the probability 
density of the wave packet is about (Ak) 2d . 
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